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Summary  of  accomplishments 

Energy  transfer  reactions  that  could  potentially  limit  scaling  of  the  electric  oxygen 
iodine  laser  (EOIL)  were  examined.  Quenching  of  excited  iodine  atoms  (I*)  by  atomic 
oxygen  has  been  identified  as  a  significant  energy  loss  channel.  The  rate  constant  for  this 
process  was  characterized  over  the  temperature  range  from  295  to  360  K. 

Quenching  of  singlet  oxygen  (02(a))  in  moderate  pressure  discharges  is  a  process 
that  exhibits  a  non-linear  pressure  dependence.  The  reactions  responsible  for  this 
behavior  have  been  probed  by  examining  the  deactivation  of  CL(a)  in  the  presence  of 
Oo(X)  and  O  atoms.  Rapid  quenching  was  observed  when  both  Oi(X)  and  O  were 
present,  suggestive  of  a  three-body  process.  However,  a  detailed  kinetic  analysis 
indicates  that  a  reaction  product  (possibly  vibrationally  excited  ozone)  is  the  primary 
quenching  agent. 

Detailed  models  of  the  post-discharge  kinetics  of  EOIL  were  fonnulated  and 
explored.  The  relative  importance  for  various  energy  loss  processes  was  evaluated,  and 
weaknesses  in  the  existing  rate  constant  database  were  identified.  An  extended  kinetic 
model  of  the  iodine  dissociation  process  was  developed.  This  model  is  applicable  to  both 
electric  and  chemical  oxygen  iodine  lasers. 


Publications  resulting  from  this  research  program: 

V.  N.  Azyazov,  P.  A.  Mikheyev,  D.  J.  Postell  and  M.  C.  Heaven,  Chem.  Phys.  Lett. 
Accepted,  Sept  2009 

"02(a  A)  quenching  in  the  O/O2/O3  system  ” 

P.  A.  Mikheyev,  D.  J.  Postell  and  M.  C.  Heaven,  J.  App.  Phys.  105  (2009)  094911 
“Temperature  dependence  of  the  O  +  I(  P 1/2)  —>  O  +  I(  P 3/2)  quenching  rate  constant ” 

V.  N.  Azyazov,  S.  Yu.  Pichugin  and  M.  C.  Heaven,  J.  Chem.  Phys.  130  (2009)  104306 
“On  the  dissociation  off  by  02<a!  A):  Pathways  involving  the  excited  species  f(A  '3ll2W 
Adi  1  „),  f(X1X,v)  and  02(aIA,v)” 

V.  N.  Azyazov,  S.  Yu.  Pichugin  and  M.  C.  Heaven,  Proc.  SPIE,  vol.  7131,  paper  22 
(2008) 

“Multi-pathway  f  dissociation  model  for  COIL  ” 
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1.  Temperature  dependence  of  the  O  +  I(2Pi/2)  —  o  +  I(2P3/2)  quenching  rate 
constant 

Discharge  singlet  oxygen  generators  add  complexity  to  the  chemistry  of  the 
oxygen  iodine  laser  as  oxygen  atoms  are  not  produced  in  the  chemically  driven  system. 
The  oxygen  atoms  in  the  post  discharge  flow  were  shown  to  effect  the  performance  of  the 
laser  to  such  a  great  extent  that  removal  of  the  oxygen  was  required  in  order  to  achieve 
gain  on  the  I*-I  transition  [1-4].  Removal  was  accomplished  by  addition  of  N02  to  the 
flow.  Several  previous  studies  have  addressed  the  collisional  deactivation  of  I*  by  ground 
state  oxygen  atoms  [1-4]. 

I*(2P1/2)  +  0(3P)  -  I(2P3/2)  +  0(3P) 

Kinetic  measurements  that  were  designed  to  isolate  the  quenching  reaction  were 
conducted  at  room  temperature  [5,6].  An  additional  complication  that  arises  with  use  of 
the  discharge  is  that  the  excess  energy  goes  into  heating  of  the  post-discharge  flow. 
Experiments  and  modeling  of  the  discharge  driven  laser  indicate  temperatures  well  above 
350  K  for  most  conditions.  Prior  to  the  work  reported  here,  no  systematic  studies  of  the 
temperature  dependence  of  this  rate  coefficient  had  been  made.  The  value  given  by  the 
recent  room  temperature  studies  seems  to  aid  in  modeling  of  EOIL,  but  measurements  of 
this  rate  coefficient  at  elevated  temperatures  were  needed.  In  addition,  temperature 
dependent  studies  also  help  determine  the  validity  of  the  previous  hypothesis  that  the 
quenching  process  occurs  via  a  curve  crossing  mechanism. 

The  rate  constant  for  quenching  of  I(2Pi/2)  by  0(3P)  was  been  measured  for  the 
temperature  range  295-360K.  Pulsed  laser  photolysis  of  mixtures  N20  and  I2  was  used  to 
examine  the  kinetics.  0(3P)  atoms  were  produced  by  the  photo-initiated  reaction 
sequence  N20+hv  — >  0(]D)+N2,  0('D)+M  —>  0(3P)+M,  while  singlet  oxygen  was 
generated  by  the  secondary  reaction  0(1D)+N20  — >  02(a' A)+N2.  Iodine  atoms  were 
produced  by  I2  photodissociation  and  from  the  secondary  reactions  of  I2  with  0(  P) 
atoms.  Subsequent  excitation  of  I  by  02(a1A)  led  to  I(2Pi/2)  formation,  with  I(2Pi/2) 
concentrations  monitored  using  time -resolved  1315  nm  emission.  The  rate  constant  of  the 
quenching  process  was  determined  by  fitting  a  kinetic  model  to  the  observed  emission 
traces.  Special  attention  was  given  to  the  gas  flow  conditions. 

A  modest  temperature  dependence  of  the  quenching  rate  constant  was  found,  and 
this  was  well  represented  by  the  equation  k  =  (6.5±1.0)xl0'12(T/300)1 76  cmV1.  The 
room  temperature  rate  constant  differs  from  previous  experimental  determinations.  The 
temperature  dependent  rate  constant  agrees  well  with  the  value  obtained  by 
computational  modeling  of  the  data  from  discharge-driven  oxygen  iodine  lasers  at  a 
nominal  temperature  of  375  K. 

2.  Quenching  of  02(a)  in  the  presence  of  O  and  02 

The  maximum  02(a’A)  yields  in  discharge  systems  obtained  so  far  (15-20  %) 
have  been  for  relatively  low  oxygen  pressures.  In  this  case  the  concentration  of  the  active 
components  in  the  flow  is  small,  resulting  in  a  low  gain  coefficient.  To  produce  a  high- 
power  oxygen  iodine  laser  system  it  will  be  necessary  to  significantly  increase  the 
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operating  pressure  of  the  device.  Braginskiy  et  al.  [7]  observed  that  the  Oita1  A)  yield  at 
the  exit  of  discharge  singlet  oxygen  generator  decreased  with  increasing  oxygen  pressure. 
They  suggested  that  faster  deactivation  of  singlet  oxygen  was  encountered  at  higher 
pressures  due  to  the  three -body  quenching  process. 

0(3P)  +  02(alA)  +  02  ->  0(3P)  +  02  +  02 

If  this  is  correct,  the  pressure  scaling  of  CPA/1  A)  densities  will  be  limited  by  three -body 
quenching  if  the  production  of  O  atoms  is  not  adequately  controlled.  To  explain  the 
singlet  oxygen  dynamics  in  the  post-discharge  flow,  Braginskiy  et  al.  [7]  proposed  that 
the  rate  constant  for  the  three-body  deactivation  was  in  the  range  k2 b  =  (l-3)xl0  32  cm6/s. 

In  the  present  work,  248  nm  laser  photolysis  of  Os/CF/Ar/CCF  mixtures  was 
employed  to  study  the  kinetics  of  CBla’A)  in  the  presence  of  oxygen  atoms,  O2  and 
ozone.  UV  photolysis  of  ozone  yields  the  products: 

O3  +  hv  ->  O^D)  +  02{aA) 

->  0(3P)  +  02(X3E), 

where  the  singlet  oxygen  yield  is  90  %.  The  02(a{A)  removal  rate  was  monitored  by 
observing  the  time  evolution  of  the  O2  emission  at  a  wavelength  of  A.=1268  nm.  The 
temporal  profiles  of  the  oxygen  atom  concentrations  were  monitored  by  means  of  the 
O+NO+M  chemiluminescent  reaction. 

Fast  and  slow  decays  of  Oita' A)  were  observed.  The  fast  decay  was  observed  when 
O  atoms  were  present  in  the  system.  An  attempt  was  made  to  model  the  fast  decay  using 
the  three  body  quenching  process  (^(a’A)  +  O  +  O2  — >  202  +  O  suggested  by  Braginskiy 

o  1 

et  al.  [7]  The  near  gas  kinetic  rate  constant  obtained  from  this  analysis  ((l.l±0.1)xl0~ 
cm6/s)  was  inconsistent  with  data  obtained  from  flowing  afterglow  experiments, 
indicating  that  additional  quenching  species  are  generated  by  the  ozone  photochemistry. 
Rate  constants  for  quenching  of  CFl^'Z)  by  CO2  ((6.1±0.5)xl0  13  cm3/s)  and  O3 
((1.9±0.2)xl0"  cm  /s)  were  measured  as  a  test  of  the  kinetic  analysis  techniques. 
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3.  Post  discharge  kinetics  in  He/02/N0/l2  mixtures 

The  relative  importance  of  various  known  reactions  in  the  flowing  afterglow  of 
electrically  discharged  He/02/N0/I2  mixtures  was  evaluated  for  the  experimental 
conditions  and  the  model  described  by  Palla  et  al.  [1].  The  results  and  modeling  of 
references  were  also  taken  into  account.  Other  rate  constants  were  taken  from 
references [8- 12].  A  kinetic  model  based  on  most  of  the  reactions  that  could  be  found  in 
literature  was  developed  and  rates  of  the  reactions  were  compared  straightforwardly  to 
detennine  which  are  most  important.  The  gas  temperature  after  the  discharge  was 
assumed  to  be  415  °K  as  that  value  provided  a  good  agreement  with  [O]  concentrations 
obtained  by  modeling  in  reference  [1].  The  transport  time  from  the  discharge  exit  to  the 
point  where  N02  was  injected  was  derived  from  the  data  in  [1]  and  assumed  to  be  30  ms. 
The  temperature  at  this  point  was  assumed  to  be  375  °K.  Then,  after  a  further  10  ms  of 
transport  time,  I2  was  injected  at  360  °K.  The  reactions  were  monitored  for  another 
10  ms,  until  I2  dissociation  was  complete.  In  the  EOIL  device  modeled,  Ar  or  N2  was 
injected,  prior  to  a  supersonic  expansion  [1].  After  that  point  the  concentrations  of  02  and 
I*  were  reduced  by  a  factor  of  about  three,  considerably  reducing  the  quenching  rates. 

Reactions  considered  in  the  model  are  listed  in  table  1 . 


Table  1. 


k 

Reaction 

Rate  constant  (cm3(cm6)s‘1) 

Ref 

1 

02(‘A) 

+ 

O^A) 

- 

OA1!) 

+ 

02 

9.8xl0’28TJ'8exp(700/T) 
9x10  17exp(-560/T) 
2.7xl0~17 

[1,5,6] 

[13] 

[14] 

2 

02(‘A) 

+ 

O^A) 

- 

o2 

+ 

02 

1.7xl0'lv 

9xl0“17exp(-560/T) 

[U3] 

3 

02(‘A) 

+ 

NO 

NO 

+ 

02 

8.5xl0‘17 

[15] 

4 

02(‘A) 

+ 

o2+o 

- 

o2+  02 

+ 

0 

1.0x1  O'32 

2.5xl0~32 

[1,15, 

16] 

5 

02(‘A) 

+ 

02 

- 

02 

+ 

02 

8.2xl0‘iy 

(1.58-2.2)xl0"18(T/300)°'8 

3  x  1 0" 1 8exp(- 200/T) 

[1] 

[7] 

6 

02(‘A) 

+ 

0 

- 

02 

+ 

0 

7xl0‘lb 

2xl0'16 

[1,7] 

7 

20 

+ 

0 

o2(‘A) 

+ 

0 

6.93xl0~35(T/300)'063 

[13] 

8 

20 

+ 

02 

- 

02(iA) 

+ 

0 

1.93xl0'“(T/300)"Ub3 

3 . 8  x  1 0"3 1  xT" 1  exp(- 1 7  0/T) 

[U3] 

9 

20 

+ 

o/A) 

-> 

o/A) 

+ 

02 

4.5xl0‘34exp(630/T) 

[1] 

10 

20 

+ 

He 

-> 

o2(1a) 

+ 

He 

9.88xl0‘35 

[13] 

11 

20 

+ 

0 

- 

02 

+ 

0 

9.21x1 0 ‘34(T/300)"ut’3 
4.5xl0‘34exp(630/T) 

[1,12, 

13] 

5 


12 

20 

+ 

02 

- 

02 

+ 

02 

2.56xl0‘24(T7300)'uw 

4.5xl0’34exp(630/T) 

3.34xl0'3°xT'1exp(-170/T) 

[1,7, 

13] 

13 

20 

+ 

He 

- 

02 

+ 

He 

lx  10" 

4.5xl0’34exp(630/T) 

[U3] 

14 

0 

+ 

02  +He 

- 

03 

+ 

He 

5.1xlO'27T"2'8 

3.4xl0'34(T/300)’12 

[U3] 

15 

0 

+ 

202 

- 

03 

+ 

o2 

5.1xlO‘2/T’24i 

6.2xl0’34(T/300)'2' 

[1,7] 

16 

20 

+ 

o2 

- 

03 

+ 

0 

4.5x10’ J4exp(630/T) 
2.5x1 0’34exp(345/T) 

[1,7] 

17 

0 

+ 

02  +  02(1A) 

02(lA) 

+ 

o3 

5.  lx  10'"  Y2’8 

[1] 

18 

o3 

+ 

(hi1  A) 

-»■ 

202 

+ 

0 

5.2x10  nexp(-2840/T) 

[12] 

19 

o3 

+ 

0 

- 

02(lA) 

+ 

o2 

2.4xl0’13exp(-2060/T) 

2.0xl0’nexp(-2280/T) 

1 .0x10" 1 1exp(-2300/T) 

[7] 

[8] 

20 

o3 

+ 

0 

02 

+ 

o2 

8xl0‘12exp(-2060/T) 

[7,13] 

21 

o3 

+ 

0 

Ool1!) 

+ 

0 

8xl0’14exp(-2060/T) 

[15] 

22 

0 

+ 

no+o2 

-»■ 

N02 

+ 

02 

4.68x10’2ST’13 

[1] 

23 

0 

+ 

NO+He 

no2 

+ 

He 

2.08x1 0'28T'I2> 

[1] 

24 

0 

+ 

no2 

- 

NO 

+ 

o2 

6 . 5  x  1 0’ 1 2exp(  1 20/T) 
9.7xl0’12 

[1]  [15] 

25 

NO 

+ 

o3 

N02 

+ 

o2 

1 .4x  1 0’ 1 2exp(- 131 0/T) 

[12] 

26 

02(12) 

+ 

0 

- 

o2(lA) 

+ 

0 

7.2xl0’14 

7.2x10’14x(T/300)°-5 

[12,13] 

27 

02(12) 

+ 

0 

- 

02 

+ 

0 

8x10° 

8x10’15x(T/300)°-5 

1.5xl0'16 

[1]  [13] 
[7] 

28 

02(12) 

+ 

o3 

- 

202 

+ 

0 

7.33xl0'12(T/300)U2’ 

1.5xl0‘u 

[1,7, 

13] [8] 

29 

02(1S) 

+ 

o3 

- 

o2(1A) 

+ 

o3 

7.33x10’12(T/300)u'5 

7.1xl0‘12 

3.3xl0‘12 

[1,7, 

13] 

30 

02(1S) 

+ 

03 

- 

02 

+ 

o3 

7.33x10’12(T/300)U3 

3.3xl0‘12 

[7,13] 

[1] 

31 

02(12) 

+ 

02 

- 

o2(lA) 

+ 

o2 

3.6xl0’lv(T/300)a5 

3.710’17 

[U3] 

32 

o2(1s) 

+ 

02 

- 

02 

+ 

o2 

4x10'18x(T/300)0-5 

1.5xl0'16 

4.3x  10’22xT’2'4exp(-24 1/T) 

[7,13] 

[11] 

33 

20 

+ 

02 

OA1!) 

+ 

02 

7.6x1  0’32xT"  'exp(- 1 7  0/T) 

[7] 

34 

I2 

+ 

0 

-»■ 

10 

+ 

I 

1.25xl0"lu 

[12] 

35 

10 

+ 

0 

02 

+ 

I 

1.4xlO'10 

[12] 

36 

10 

+ 

10 

02 

+ 

I2 

5.4x1 0"cxp(180/T) 

[12] 

6 


I  + 


38  I* 

39  _ 

40  I* 


I _ + 

I  +" 


44  I* 


46  I* 

47  h*  +" 

48  I2*  + 

49  O^Z)  + 

50~  O^Z)  T 

5^_0/2)_^ 

52__02(^L^ 

53  02(1A)  + 

54  O  + 
55~  02('X)  T 


Q2(1A) 

02 

NO 


NO+M 

INO 

02(1A) 

02('A) 


o^a) 

02 


I* 

02(1A) 

NO 


+  02 

-+~r 


o2Cz) 

02('A) 


O^A) 

O^A) 

o^a) 

02 


+  M 
+  h 


+  21 
+  I2 

+  21 

±JT 

+  NO 
+  He 


2.3x10  /T 

3.1x1 0"8xT‘'exp(-403/T) 
1.2x10'' 

3.8x10" 

1 .  4x  1  O' 1 3exp(  1 600/T) 
1.7x10"' 
1.6xl014 
1.8x10"' 
2.6x10"° 
l.lxl  0  1 ' 

4x  1 0~24xT"3  8xexp(7 00/T) 
1.1x10"' 

6.5xl0"-x(T/300)1Vb 

_ 8xl0~12 _ 

3x10"° 

4.9x10" 2 
3.8x10"' 

_ 4xl0~12 _ 

2.3xl0‘n 

6.0xl0"4 

_ l.OxlO'17 _ 

3x10'5x(v)/2R 

4x10'4x(v)/2R 

(1-3)x10'3x(v)/2R 

1x10'2x(v)/2R 


o2 a,  no  and  o3p  at  4 15  K 


The  kinetics  for  different 
species  were  analyzed  separately  to 
understand  the  importance  of 
reaction  rates  that  cause  their 
production  and  removal  during 
specified  time  intervals.  02(a'A), 
O^b'l),  03,  O,  I*,  I,  NO,  N 62  and 
O2  are  denoted  in  the  figures  as  o2a, 
o2b,  o3,  o3p,  iex,  ix,  no,  no2  and 
o2x  . 

Initial  values  for  the 
concentrations  right  after  the 
discharge  were  chosen  to  be  close  to 
those  given  in  [1].  Specifically, 
[O2]=4.8x1016,  [O2(1A)]=6.9x1015, 

[O2(1E)]=6.9x1014,  [O3]=0,  [0]= 

4.8xl015,  and  [NO]=2.9xl015  cm'3 
at  415  °K. 

The  evolution  of  [O2CA)],  [O]  and  [NO]  during  the  30  ms  transport  is  represented 
in  Figure  1.  The  number  densities  for  this  group  are  of  the  same  order  of  magnitude 
(10  cm'  )  during  first  30  ms  of  transportation.  NO  added  to  the  discharge  favors  singlet 
oxygen  production  [1,15]  and  provides  the  fastest  channel  of  O  atom  removal  in  the 
reaction  sequence  22-24.  O3,  N02  and  02('Z)  are  also  present  in  the  gas  stream  from  the 
dischsrge.  These  species  have  number  densities  two  orders  of  magnitude  lower  than  those 
of  Fig.  1.  The  evolution  for  this  second  group  during  transport  is  shown  in  Fig.  2. 

o2b,  o3  and  no2  at  415  K 

Singlet  delta  oxygen. 

The  rates  for  quenching  of  02( 1  A) 
are  represented  in  Fig.  3.  The  reactions 
are  listed  in  the  caption  according  to  their 
rates  in  descending  order.  They  are 
reactions  4,  18,  6,  1,  3,  53,  2,  5  in  Table 
1.  Rates  2  and  5  contribute  just  13  %  of 
the  quenching,  and  only  at  the  end  of 
30  ms  time  interval.  Hence,  these 
reactions  may  be  neglected  at  shorter 
times.  Although  [O3]  is  quite  low  its 
contribution  to  02(1A)  quenching  is  one 
of  the  largest.  It  should  be  noted  that  the 
rate  constants  of  the  most  important 
reactions  4  and  6  are  known  poorly  and 
different  authors  use  quite  different 
values,  as  seen  from  the  table. 


Figure  2.  Evolution  of  02(b),  N02 
and  O3  during  30  ms  transport. 


Time 

Figure  1.  Evolution  of  02(a),  NO  and 

•5 

0(  P)  during  30  ms  transport. 
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Wall  recombination  was 
considered  in  refs  [15,17].  The  rate 
coefficient  for  this  reaction  is 
estimated  by  kw=y(v)/2R,  where  10‘6 
7=3x10°  [15]  is  the  recombination 
probability,  (v)«3xl04  cm/s  is  the 
mean  oxygen  atom  thermal  g 
velocity,  and  2R=5  cm  is  the  tube  J 
diameter  [1].  As  can  be  seen  in  Fig.  i 
3,  wall  recombination  should  be 
taken  into  account  on  this  time  scale 
because  at  the  beginning  its 

io*5 

contribution  is  5%  and  at  the  end  it 
has  increased  to  20%. 

Some  Ooi’A)  is  produced  in 

i 

the  afterglow,  mostly  due  to  C>2(  2) 
conversion,  resulting  in  a  small  rise 
of  [(At1  A)]  right  after  discharge. 

Rates  for  the  reactions  7,  8,  and  31  are  smaller  than  3X1014  cm'V1  and  may  be  safely 
excluded  from  the  model  as  their  contribution  to  CF^A)  production  is  less  than  3%. 
Different  authors  [7,8,18]  use  rate  coefficients  for  reaction  19  with  2  orders  of  magnitude 
difference.  If  the  largest  estimate  is  correct,  this  reaction  should  be  taken  into  account. 

For  the  experimental  conditions  of  [1]  the  overall  loss  of  Ch( 1  A)  during  first 
30  ms  after  discharge  is  about  7%.  This  figure  remains  the  same  regardless  of  the 
presence  or  absence  of  NO  in  the  system.  The  reason  is  that  in  the  absence  of  NO,  despite 
a  larger  quenching  by  O,  some  Oo^A)  is  produced  by  reaction  10.  However,  for  the 
higher  yield  of  [Oo^A)]  obtained  in  [15]  and  with  O  atoms  removed  by  HgO  in  the 
discharge,  the  presence  of  NO  facilitates  02^  A)  quenching. 

Oxygen  atoms 

In  the  presence  of  NO  the  fastest  rates  of  O  atoms  removal  are  from  reactions  24, 
23,  22,  53,  13,  14,  12,  and  15.  The  reaction  rates  are  in  the  range  1 0 1 5- 1 0 1 7  cm'V1.  It  is 
important  to  note  that  the  next  largest  removal  rate  is  provided  by  wall  recombination 
with  7=4x1  O'4  [17]  which  must  be  accounted  for  in  the  modeling.  In  the  absence  of  NO 
the  rate  of  heterogeneous  recombination  dominates.  Reactions  18,  20,  16,  and  10,  with 
rates  are  in  the  range  (0. 1  -5)x  1 0 15  cm'V1  are  also  important  for  modeling,  as  they 
determine  [O]  in  the  absence  of  NO. 

The  rates  of  reactions  7-9,  11,  19,  21  are  less  than  2x10 14  cm'V1  and  their  sum  is 
less  than  5X1014  cm'3 s'1.  Consequently,  they  may  be  excluded  from  the  model  (bearing  in 
mind  the  problem  with  the  rate  constant  for  reaction  19)  as  their  contribution  to  the  total 
[O]  removal  is  about  1%. 


o2 a  quenching  at  415  K 


3.  02(a)  quenching  rates. 
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o2b  quenching  at 415  K 


O2C Z)  and  O3 

In  the  absence  of 
iodine,  CF^E)  is  produced  in 
reactions  1  (pooling)  and  33. 

The  latter  contributes  less  than 
5%  and  may  be  safely 
discarded.  Removal  of  CFfS) 
occurs  via  reactions  26,  51,  28, 

55,  27,  29,  30  corresponding  to 
their  descending  rates  and  the 
order  they  appear  in  Fig.  4. 

Reactions  31,  32,  52  may  be 
neglected,  being  two  orders  of 
magnitude  slower.  After  the 
removal  of  most  of  the  O 
atoms,  quenching  of  CFt'l)  by 
NO  and  wall  collisions 
(Y=lxl0‘2)  [18]  are  the 

dominant  removal  pathways. 

The  concentration  of  ozone  generated  by  the  discharge  is  small  [1,13].  Ozone 
production  in  the  afterglow  occurs  via  reactions  14-17  and  it  is  removed  by  reactions  25, 
18,  19,  and  28.  Reactions  17  and  19  contribute  less  then  3%  to  the  rate  and  may  be 
discarded.  After  a  brief  increase,  [O3]  decreases  during  transport,  and  after  30  ms  drops  to 
~5  x  1 012  cm'3  (Fig.  2). 


Behavior  of  the  system  after  NO2  and  I2  injection 

After  NO2  injection  the  O 
atoms  are  rapidly  eliminated  from 
the  system  and  slow  02(1A) 
quenching  occurs  due  to  reactions 
3  (NO  quenching),  1  (pooling),  53 
(wall),  2  (self  quenching)  and  5 
(by  O2)  resulting  in  less  than  2% 
decrease  in  [CF^A)]  during  10  ms 
before  F  is  injected. 

After  F  injection,  the  loss 
of  02(1A)  substantially  increases 
due  to  I*  quenching  reactions. 

Fig.  5  shows  that  F  dissociation  is 
complete  after  10  ms.  Details  of 
the  F  dissociation  mechanism  are 
not  considered  here  and  the 
process  is  treated  empirically 
using  reactions  40  and  47. 

Experiments  and  modeling 
[1]  have  shown  that  lasing  can  be 


Concentrations  after  12  injection  at  360  K 


Figure  5.  Concentrations  after  F  injection. 
[02(1A)]  is  scaled  to  fit  the  figure. 
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I  *  loss  rates  at  360  K 


achieved  only  when  most  of 
the  O  atoms  are  removed  from 
the  system.  The  present 
calculations  indicate  that  the 
dissociation  of  F  proceeds 
without  significant 

contributions  from  reactions 
34  and  35.  Even  at  a  low 
level,  the  presence  of  O  atoms 
decreases  [02(1A)]  during  E 
dissociation. 

Rates  of  I*  loss  are 
depicted  in  Fig.  6.  Reactions 
I*+wall,  44,  45,  39,  40 

correspond  to  the  rates  in 
decreasing  order.  As  there  is  a 
near  equilibrium  between 
02(1A),  02,  I  and  I*, 

irreversible  loss  of  I*  is 
equivalent  to  loss  of  02(1A).  This  makes  wall  quenching  of  I*  an  important  process.  This 
situation  was  considered  in  detail  in  reference  [19],  and  it  was  found  that,  for  the  typical 
conditions  of  a  chemical  oxygen-iodine  laser,  wall  loss  may  be  assessed  by  assigning  a 
wall  loss  probability  for  Ch^A)  of  the  order  [I]/[C>2],  assuming  that  every  I*+wall 
collision  results  in  deactivation.  With  the  increase  of  [I]  during  dissociation  this  value 
changes  and  under  the  conditions  of  reference  [1]  the  mean  estimated  wall  loss 
probability  is  2x10°.  Fig.  6  shows  this  to  be  the  dominant  loss  process  after  F  addition. 

The  rate  constant  of  reaction  41  used  in  [1]  is  an  order  of  magnitude  larger  than  in 
[14],  but  its  rate  is,  nevertheless  low  and  it  may  be  disregarded.  The  rate  of  NO  assisted  I 
recombination  due  to  reactions  42,  43  is  smaller  than  rate  of  reaction  41  and  may  also  be 
neglected.  In  the  absence  of  O  atoms  reaction  46  does  not  contribute  to  I*  quenching. 

Conclusions 

Of  the  reactions  listed  in  Table  1,  15  have  been  identified  that  are  not  of 
importance  in  modeling  EOIF  systems.  These  are  reactions  7-9,  1 1,  17,  19,  21,  31-33,  41- 
43,  52.  The  rate  coefficient  of  reaction  19  differs  by  two  orders  of  magnitude  in  the 
literature.  If  the  true  value  is  close  to  the  upper  limit,  reaction  19  should  remain  in  the 
model.  Heterogeneous  quenching  is  identified  as  being  important  in  the  post  discharge 
kinetics. 


Time 

Figure  6.  Foss  rates  after  F  injection. 
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4.  On  the  dissociation  of  I2  by  O^a'A):  Pathways  involving  the  excited  species 
I2(A'3n2u,  A3nlu),  I2(X!2:,i))  and  02(a*A,  n) 

Kinetic  studies  were  carried  out  to  explore  the  role  of  the  excited  species  I2(A'  /72u, 
A377iu),  I2(X' A,  o)  and  Oita1  A,  o)  in  the  dissociation  of  I2  by  singlet  oxygen.  A  flow  tube 
apparatus  that  utilized  a  chemical  singlet  oxygen  generator  was  used  to  measure  the  I2 
dissociation  rate  in  O £)l\i  mixtures.  Vibrationally  excited  I2(X)  is  thought  to  be  a 
significant  intermediate  in  the  dissociation  process.  Excitation  probabilities  (yv)  for 
population  of  the  u-th  I2(X)  vibrational  level  in  the  reaction  hiX)+\(  P\a)  —*■  I2(X 
d>10)+I(2P3/2)  were  estimated  based  on  a  comparison  of  calculated  populations  with 
experimentally  determined  values.  Satisfactory  agreement  with  the  experimental  data 
was  achieved  for  total  excitation  probabilities  partitioned  in  two  ranges,  such  that  /25<y<#7 

47  23 

=  X  Yo  ~  0-1  and  7~/5<o<25  =  Z  ~  0-9.  A  multi-pathway  I2  dissociation  model  was 

u=25  u=15 

developed  in  which  the  intermediates  are  I2(A,3772u,  A  'fl\ u)  and  h(Xjj).  It  was  shown  that 
the  iodine  dissociation  process  passes  predominantly  through  the  I2(A,3/72u,  A3I7lu) 
intermediate.  These  states  are  populated  mainly  by  collisions  of  I2  with  vibrationally 
excited  O 2(aA,  v )  when  the  mole  fraction  of  I2  is  small  (rjj  <  1%),  and  by  collisions  of 
I2(X15</j<23)  with  O i{a  A)  for  higher  concentrations  ( i]j  >\%).  The  revised  kinetic 
model  is  defined  below  in  Table  4.1 
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Table  4.1.  The  extended  kinetics  model  for  COIL 


#  Reaction 


1 

2 

3 

4 

5 

6 

7 

8 
9 

14 

15 

16 

17 

18 
21 
22 

24 

25 

26 

27 

28 

29.1 

29.2 

29.3 

29.4 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


Oo  (a)  +  O2O2)  — > 
02(a)  +  02(a)  — > 
02(a)  +  02(a)  —> 
O  2(b)  +  02(X)- 
02{b)  +  U20  -h 
O  2(b)  +  Cl2  - 
O 2(b)  +  H2O2  - 
O  2(b)  +  He  -> 
O  2(b)  +  C02  - 

02(a)  +  02(X)^ 
O  2(a)  +  H20  — s 
02(a)  +  Cl2  -h 
02(a)  +  H2O2  — ■ 
02(a)  +  He  — » 
h(X)  +  O 2(b)  -» 
I2W  +  02(b)  -> 
h(X)  +  O 2(b)  -» 
HA,  A')  +  02(a) 


O  2(b)  +  02(X) 
02(a)  +  02(X) 
02(X)+  02(X) 
02(a)  +  02(X) 
02(a,  0 )  +  H20 
02(a)  +  Cl2 
•  02(a)  +  H202 
02(a)  +  He 
->  02(a,  0 )  +  C02 
o2(X)  +  02(X) 

02(X)  +  h2o 

02(X)  +  Cl2 
02(X)  +  H2O2 
02(X)  +  He 
I  +  I  +  02(X) 
I2W  +  02(a) 
HA, A')  +  02(X) 
->  I  +  I  +  02(X) 
-  UB)+o2(X) 


I 


HA, A')  +  02(a) 

I2(5)  +  02  — *  I  +  I  +  02 
I2(A)  +  O2(X)->l2(20+O2(a) 

HA  )  +  H20  -»•  I2  +  H20 
ha  )  +  No  21  +  N2 
h(A  )  +  C02  ->  I2  +  co2 
HA  ')  +  He  -*■  21  +  He 
1 2(5)  —  I  +  I 
HE)  HX)  +  hv 
I2(X)  +  02(a)  HX,v)  +  02(X) 

I2(X)  +  I*  -»■  HX,v  >10)  + 
b(X,v>25)  +  02(a)  ->  I  +  I  +  02(X) 

I2(X,v)  +  02^  I2(X,o-l)  +  02 
HX,v)  +  H20  -»■  I2(X,  o-l)+  H20 
•  l2(Xv-l)  +  He 
I2(X,o-l)  +  N2 
►  I2(X,o-l)  +  C02 
I*  +  02(X) 

I 
I 
I 
I 
I 


HX,v)  +  He 
HX,v)  +  N2 
HX,v)  +  C02 


I 

I* 

I 

I* 

I* 

I* 


+ 

+ 

+ 

+ 

+ 

+ 


02(a) 
02(X) 
o  2(a) 
02(X) 
02(a) 
02(a) 


+  02(a) 
+  02(X) 
+  02(X) 
+  O  2(b) 
+  02(a) 


Rate  constant  (7=300  K), 

_ cm3/s  _ 

2.5xl0"lv 

0 

1.7xl0"17 

3.9xl0"17 

6.7xl0"12 

2.0xl015 

3.3xl0"13 

l.OxlO"17 

6xl0"13 

1.6xl0~18 

4.0xl0"18 

6.0xl0~18 

0 

8.0x10 21 
3.3x10" 
2.5x10" 

0 

3x10  11 
<10‘12 
6x10  11 
6.3xl0~12 
3.4xl0~12 
3.5xl0~14 
8.5xl0"13 
9.4xl0~15 
6xl05 s" 

6x10s s" 
<5.0xl0"16 
3.8x10" 

3x10  11 
ox2.7xl0  12 
ox6xl0  12 
ox3.9xl012 
ox3.4xl012 
ox3.4xl0"2 
7.8x10" 
2.7x10" 
l.OxlO"15 
0 

l.lxlO"13 

l.lxlO'13 
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46 

47 

48 

49 

50 

51 

53 

54 

56 

57 

58 

59 

60 
61 
62 

70 

71 

72 

73 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

89 

90 

91 

92 

93 

94 

95 

96 

97 
101 
102 
103 


I* 

I* 

I* 

I* 

I* 

I* 

I* 

I* 

I* 

h 


+ 

+ 

+ 

+ 

+ 

+ 


O  2(a) 

I 

H20 

h2o2 

He 

02(X) 


->  I  +  O2(30 

->  I  +  I 

-»•  I  +  h2o 

-»•  I  +  h2o2 

— >  I  +  He 

-»•  I  +  02(X) 

—>  I  +  hv 
IC1  +  Cl 
-»•  I2  +  Cl 
— >  IC1  +  I 

—  ci2  +  I 

->  I2(X)  +  I2(3Q 

—  I2(B)  +  b(X) 
->  I2(X)  +  He 

—  I2®  +  02(X) 
02(o-l)  +  02 
02(o-l)  +  H20 
02(o-l)  +  He 
02(u-l)  +  N2 


+  Cl2 
+  IC1 
+  Cl 
IC1  +  Cl 
I  +  I  +  l2(X) 

I*  +  I  +  l2(X) 

I  +  I  +  He 

I  +  I  +  O2(30 

02(w)  +  02  -» 

O  2(v)  +  H20  — > 

02(o)  +  He  — > 

02(o)  +  N2  -> 

02(X,2)  +  O2(X,0) 

02(X,3)  +  O2(X,0) 

02(X4)  +  O2(X,0) 

O  2(ci,  1)  +  C02  — > 

02(a,2)  +  C02  — > 

02(a,3)  +  C02  — > 

02(&,l)  +  C02-> 

02(b,l)  +  C02  — > 

02{b,3)  +  C02  — > 

02(X1)  +  H20(000) 

02(a,  1)  +  O2(36,0)  - 
02(a,2)  +  O2(X,0)  - 
O2(b,\)  +  O2(X,0)  - 
02(b,l)  +  O2(36,0)  - 
02{b,3)  +  O2(36,0)  - 
H2O(010)  +  H20  -> 

02(a,  1)  +  I2(3Q  — » 

02(a,2)  +  I2(3Q  — > 

02(a,3)  +  1 2  (36)  — > 

02(o)  +  I2(36,  1 1<  o  <24) 

02(a,  1)  +  I2(36,  o  >15) 

02(o,2)  +  I2(X,d>8)  ->  02(X)+2I 


->  02(36,1)  +  02(X,l)a 
->  02(X,2)  +  02(36,1) 

->  02(X3)  +  02(X1) 
02(0,0)  +  C02(o) 

02(a,  1)  +  C02(o) 

02(o,2)  +  C02(o) 

O2(b,0)  +  C02(v) 

02(b,l)  +  C02(v) 

02(b,2)  +  C02(v) 

->  O2(X,0)  +  H2O(010) 
->  02(X,1)  +  02(0,0) 

->  02(36,2)  +  02(0,0) 

->  02(36,1)  +  O2(fc,0) 

->  02(36,2)  +  O2(fc,0) 

->  02(36,3)  +  O2(fc,0) 
H20(000)  +  H20 
02(3Q  +  I2(A ') 

O2(30  + 1 2(A) 

02(X)  +  21 

►  O2(30  +  I2(A') 
02(3Q  +  21 


0 

1.6xl014 
2.0xl0~12 
2.5x10  11 
5.0xl0~18 
3.5xl0"16 
7.8  s’1 
5.5xl0~15 
1.5x10  11 
2.0xl0"1() 
8.0xl0"12 
3.6xlO’30  cm6/s 
<3.6xlO~30  cm6/s 
3.6xl0‘33  cm6/s 
3.7xl0‘32  cm6/s 
8.2xl0~19 
8.2xl0"17 
0 
0 

2.0xl0‘13 
2.6xl0"13 
2.7xl0‘13 
1.8xl0~14 
4.4xl0~14 
l.OxlO"13 
1.2xl0"12 
1.7xl0"12 
1.6xl0"12 
1.7xl0"12 
5.6xlO‘u 
3.6xlO'u 
1.52xlO"]  1 
1.7xl0’12 
1.5xlO"13 
5xlO‘n 
2xl012 
3x10" 

10'" 

1012-l0n 

1011 

10'" 


a  In  the  symbol  02(36.  v)  o  is  the  vibrational  quantum  number. 
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5.  Spectroscopic  detection  of  iodine  oxides  in  the  reactions  of  iodine  with  discharge 
excited  oxygen 

Kinetic  modeling  of  EOIL  systems  suggests  that  the  reaction  set  currently  in  use 
does  not  adequately  reproduce  the  deactivation  kinetics  of  I*.  It  appears  that  quenching 
by  some  as  yet  unidentified  species  is  involved.  We  are  exploring  the  possibility  that  IO2 
and  higher  oxides  of  iodine  may  be  responsible  for  some  of  the  quenching. 

Cavity  ring-down  spectroscopy  (CRDS)  was  used  to  observe  iodine  oxides 
formed  in  O2/I2  discharges.  Fig.  7  shows  a  typical  result  from  this  experiment.  This 
rotationally  resolved  band  has  been  identified  as  the  2-9  A  II3/2-  X~Il3/2  transition  of  10. 
The  negative  going  trace  in  Fig.  7  is  a  simulation  of  the  band  based  on  the  molecular 
constants  of  Durie  et  al.  (Can.  J.  Phys.  38,  444  (I960)),  with  a  rotational  temperature  of 
40  K.  Note  that  the  lower  vibrational  level  of  this  transition,  v”=9,  lies  5747  cm"1  above 
the  zero-point  energy.  Hence,  the  observation  of  this  transition  indicates  that  vibrational 
relaxation  in  the  jet  was  relatively  inefficient. 


In  addition  to  bands  that  could  be  assigned  to  known  states  of  10,  there  are  several  other 
bands  where  we  have  not  yet  identified  the  carrier  or  the  transition  type.  Fig.  8  shows  an 
example  of  one  of  the  new  features. 
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Figure  8.  Rotationally  resolved  spectrum  of  an  F/CF  discharge  product 

The  rotational  line  spacings  in  Fig.  8  indicate  a  relatively  large  rotational  constant, 
comparable  to  that  of  IO(X).  We  are  examining  the  possibility  that  the  bands  that  show 
structure  similar  to  that  of  Fig.  8  arise  from  the  ATT1/2-  XTI1/2  system  of  10.  The  Q=l/2 
components  of  the  A  and  X  states  have  not  been  observed  previously,  and  the  spin-orbit 
coupling  constants  are  unknown. 

In  addition  to  the  features  of  10  and  IO2,  several  new  band  systems  have  been 
observed  in  the  CRDS  traces.  An  example  of  a  low-resolution  survey  scan  is  show  in 
Fig.  9.  Further  studies  that  include  mass  spectrometric  detection  are  needed  to  identify 
the  carriers  of  the  unassigned  bands. 


0.006  J 


572  576  580 

Wavelength  /nm 


Figure  9.  Absorption  spectrum  of  the  products  from  a  discharge  through  a  mixture  of 
oxygen  and  iodine. 
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